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Abstract: We clarify the relation between the variable mT2 and the method of kine- 
matic constraints, both of which can be used for mass determination in events with 
two missing (dark matter) particles at hadron colliders. We identify a set of minimal 
kinematic constraints, including the mass shell conditions for the missing particles and 
their mother particles, as well as the constraint from the measured missing transverse 
momentum. We show that mT2 is the boundary of the mass region consistent with 
the minimal constraints. From this point of view, we also obtained a more efficient 
algorithm for calculating mT2- When more constraints are available in the events, we 
can develop more sophisticated mass determination methods starting from the 171^2 
constraint. In particular, we discuss cases when each decay chain contains two visible 
particles. 
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1. Introduction 

Many extensions beyond the standard model (SM) predict a dark matter candidate, 
that is, a massive stable particle interacting weakly with the SM particles. These 
include the minimal supersymmetric standard model (MSSM) with R-parity, universal 
extra dimension (UED) models |Q with KK-parity, little Higgs models with T-parity 
(LHT) [0] and so forth. A common feature of these models is that they all possess 
an exact parity, under which the SM particles are even and some new physical states 
are odd. Therefore the lightest parity-odd particle is stable and plays the role of the 
dark matter candidate. At a collider, the parity odd particles must be produced in 
pairs. Each of them will then go through cascade decays ending at the stable particle. 
Because it is weakly interacting, the stable particle will escape the detector without 
being detected, leaving missing energy signals. 

The Large Hadron Collider (LHC) will start collisions and collecting data very soon. 
Once large missing energy signals are detected at the LHC, it is crucial to investigate the 
properties of the new particles involved in the events, in order to identify the underlying 
theory, as well as to determine if the invisible particle is a viable dark matter candidate. 
Of particular interest are the masses of the new particles, including the invisible one. 
Due to the fact that there are always two or more invisible particles in each of such 
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events, mass determination will be a difficult task. For a hadron collider, the total 
momentum in the beam direction is not measured, making it even more challenging. 

Distributions of some simple kinematic variables such as '^Ti Ej- and M^ff 
have been used to give estimates of the masses of the new particles. However, these 
variables are mostly sensitive to the mass differences of the new particles, instead of 
the absolute mass scale. On the other hand, the total production cross section and 
the full likelihood method require knowledge of the underlying physics such as the 
matrix elements, and hence are model-dependent. Measurements of the new particle 
properties should be the ffist step towards uncovering the underlying theory instead 
of the other way around. It is therefore desirable to be able to determine masses in a 
model-independent way by using only kinematics. Traditionally, this has been done by 
looking for the edges/endpoints of various invariant mass distributions of the visible 
particles The positions of these edges/endpoints are functions of the masses of the 
particles involved in the decay chains. If the decay chains are long enough, there may be 
enough independent invariant mass endpoints involving the visible SM particles, which 
allows one to reverse the relations to obtain the masses. Nevertheless, this method 
applies to individual decay chains and can only work for the long decay chains (4 or 
more on-shell particles in a decay chain). It does not utilize all information in the events 
such as the measured missing transverse momentum. Consequently, a large number of 
events are required to distinguish and measure all edges/endpoints in order to achieve a 
reasonable determination of the masses. It is therefore important to develop new mass 
determination techniques which are more powerful and can also be applied to shorter 
decay chains for events involving invisible particles. 

Recently, along this direction, two kinds of mass determination techniques have 
been proposed. One of them utilizes the "kinematic constraints" |], 0]. Assuming 
that the event topology is known, one can try to reconstruct the kinematics event by 
event, by imposing the mass-shell constraints and the constraint from the measured 
missing transverse momentum, Depending on how many constraints that we can 
impose, the detailed methods can be different. In Ref. 0, the authors considered 
events with two identical decay chains, each containing two visible particles. Assuming 
that all intermediate particles are on-shell, the two invisible particles' momenta can be 
solved for given trial masses. Requiring the solutions to be physical, one can determine 
the masses by examining the number of solvable events for all possible trial masses. 
Another method requires longer decay chains and therefore more constraints. It is 
then possible to combine multiple events and solve directly for the masses and the 
momenta, without assuming any trial masses [0. This technique has been combined 
with the edge/endpoint method to achieve further improvement 

In a seemingly parallel approach, several authors have studied mass determination 
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Figure 1: An event with two invisible particles A^, each from a decay of a heavy particle Y. 



methods using the variable 771^2 0, which is sometimes called the stransverse mass. 
mT2 is defined event by event as a function of the invisible particle mass. Its endpoint 
or maximal value over many events, denoted by m™2^, gives an estimate of the mother 
particle's mass in the beginning of the decay chain. When the invisible particle's mass 
is unknown, one has to use a trial mass to calculate and only obtains an estimate 
of the mass difference. However, it has been shown in Ref. that if the two mother 
particles decay through three-body decays to the invisible particles, a "kink" occurs on 
the m™2^ curve as a function of the trial mass. The position of the kink is actually at the 
true value of the invisible particle mass, which allows us to simultaneously determine 
the masses of both the invisible particle and its mother particle. A generalized study 
of the kink method is available in Ref. |jll] . 

The purpose of this paper is to clarify the relation between the two mass deter- 
mination techniques, i.e., the one using kinematic constraints and the one using the 
variable mT2- An apparent difference between the two approaches is that the former 
uses the 4-momenta of the visible particles, while the latter is defined solely on the 
plane transverse to the beam direction. Nevertheless, due to the lack of total momen- 
tum measurement in the beam direction, the longitudinal momenta of the two invisible 
particles can be arbitrarily chosen, offsetting some of the information obtained from 
the visible particles' longitudinal momenta. As a consequence, 771^2 is equivalent to the 
"minimal" kinematic constraints discussed below. 

We illustrate our definition of "minimal" constraints in Fig. |l]. Two mother par- 
ticles of the same mass, my, each decays to a dark matter particle of mass rriN, plus 
some visible particles, either directly or through other on-shell particles. Since the 
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two masses are unknown, we have to assume some trial masses, denoted by hy and 
/iTV- Upstream transverse momentum (UTM) can be present but it must be known. 
The "minimal" constraints are then defined as the mass-shell conditions from /iy and 
/iTV, plus the constraint from the measured ^t- Obviously, for a given trial mass of 
the invisible particle fi^, the mother particle cannot be too light otherwise we cannot 
obtain physical momenta for the invisible particles. As we will see, we can satisfy the 
minimal constraints and obtain physical momenta if and only if fiy > ?^t2- We note 
that this fact has been implicitly used in Ref. [|l^, and we give a detailed account in this 
article. An important by-product of our discussion is that we develop a new algorithm 
to calculate mT2 which is 5-9 times as fast as the currently available program. 

Since mT2 corresponds to the minimal kinematic constraints, it can serve as a 
starting point for mass determination with more complicated topologies. The afore- 
mentioned "kink" method is an example. The m2c variable defined in Ref. [12| is 
another example, in which the authors combine mT2 with the measurement of the 
invariant mass distribution endpoint. We will give more examples in this article. 

The paper is organized as follows. In the next section, after reviewing the definition 
of mT2, we prove that it is equivalent to the minimal kinematic constraints. Deriving 
mT2 from kinematic constraints also provides us a fast algorithm for calculating mT2, 
which is presented in Section |2]^. In Section ^ we discuss some mass determination 
methods from our understanding of the relation between mjn2 and kinematic constraints. 
A more general discussion of kinematic constraints and conclusions are contained in 
Section |. 



2. mT2 from kinematic constraints 
2.1 The definition of mT2 

The definition of mT2 is motivated from the transverse mass m^, which is defined 
for events with one invisible particle at a hadron collider. In this case, the measured 
missing transverse momentum is equal to the transverse momentum of the invisible 
particle, ttit has been used, for example, in the measurement of the W boson mass in 
the decay W — ^ iu. Using the notation of the W decay, ttit is defined by 

m^ = m2 + m2 + 2(E^E^-p^-p^), (2.1) 

where mi and rrii, are respectively the masses of the lepton and the neutrino, and 
and p^ are their transverse momenta. The beams are chosen to be along the z direction, 
therefore p^ = {px,Py)- E^, are transverse energies defined by 

Et = + |p^|2, = + Ip^P. (2.2) 
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For convenience, we use a to denote the 2+1 dimensional momentum, a = {Et,Pt), 
while the 4-momentum is denoted by p. In the 2+1 dimensional notation, the transverse 
mass is given by 

ml = {ae + a,f. (2.3) 

The following relation is always satisfied for the physical momenta p^, p^ and the 
corresponding 2+1 dimensional momenta: 

{pi + Puf>{ai + a,f. (2.4) 

The equality holds if and only if £ and v have the same rapidity, which is given by 



' 2 \E-p, 

When the event contains two or more missing particles, we can no longer calculate 
the transverse mass because the transverse momentum of the individual missing particle 
is unknown. As mentioned in the Introduction, a particular interesting case is that there 
are two decay chains in the event, each ends with an invisible particle of species N . We 
further assume that each decay chain also contains a particle of species Y , decaying 
to the particle N plus some visible particles. This is illustrated in Fig. |l|, where we 
have labeled the invisible particles as 1 and 2, and summed the visible 4- momenta to 
Pa and for the two decay chains respectively. We will treat a and h as two particles 
whose masses may vary from event to event. There can be other upstream transverse 
momentum (UTM) from, for example, initial state radiation or heavier particle decays. 
However, it is important that 1 and 2 are the only invisible particles. 

Comparing with the W decay example, we see two difficulties associated with the 
above decay chains. First, m^r, the mass of the particle is a priori unknown. Second, 
only the sum of the two invisible particles' transverse momenta is measured. These 
difficulties motivated the authors of Ref. to define a quantity mT2, using a trial 
N mass (denoted by /itv) and minimizing over all possible partitions of the measured 
transverse momentum: 

m|2(Ai7v) = min [max{m^(p^, p^; /itv), m|(p|, p^; ^in)}] 

= min [max{(Q;i + a^)^, ("2 + ttfe)^}] , (2.6) 

where in the second line we have rewritten the transverse mass using the 2+1 dimen- 
sional notation. By definition, mT2 is an event-by-event quantity depending on the 
trial mass /xat. Therefore, strictly speaking it is not a variable, but a function of fiN- 
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For a given n^, we can examine the 171x2 distribution for a large number of events, 
which in general has an end point. As discussed in Ref. 0, the 771^2 end point gives 
the correct mass of the particle Y when the trial mass is equal to the true mass of the 
missing particle A^, fiN = We can therefore use mT2 to determine my if is 
known, analogous to the W mass measurement. Moreover, it has recently been shown 
rU| that, even if m^r is unknown, in some cases, when we plot the 771x2 endpoint as a 



function of the trial mass /iat, there is a kink at /iA? = mjv. Thus both m^v and my can 
be determined by studying the 771x2 distribution. 

We will discuss mass determination using 7nx2 in Section ^. Before that, we first 
give an alternative definition of mx2, using the concept of kinematic constraints. 

2.2 771x2 from minimal kinematic constraints 

By kinematic constraints, we mean two kinds of constraints imposing on the 4-momenta 
of the invisible particles: the mass shell constraints and the measured missing transverse 
momentum constraints. Specifically, for the event in Fig. |I], we can write down the 
following equations: 

2 _ 2 _ 2 

Pi — P2 ~ f^Ny 
(Pl+Pa)^ = {P2+Phf = /^y, 

Pi+P2=f,P\+Pl=f, (2.7) 

where /iy is a trial mass for the particle Y. We call this set of constraints "minimal" 
because they correspond to the shortest decay chains. Note that for a given set of 
(/iAT, /iy), the system contains only 6 equations, which are not enough for completely 
determining pi and p2- Nevertheless, Eqs. (|2.7|) still constrain the possible (/itv, /iy). 
In particular, we will shortly see that for a given fi^, Eqs. ( |2.7|) can be satisfied for 
some physical momenta pi and p2 if and only if /iy > 77ix2{fJ'N)- Here, a momentum is 
"physical" if all of its components are real and the energy component is positive. In 
other words, 7nx2{fJ'N) can be defined as the boundary of the consistent region on the 
{fiN, fJ'Y) plane, subject to the minimal constraints in Eqs. ( |2.7| ). This fact has been 
used in Ref. but without a clear proof. 



First, it is easy to show that /iy cannot go below 771x2 for a fixed fi^. For 
any (/ijv, /iy) in the consistent mass region, there exist physical pi and p2 satisfying 
Eqs. (|2.7|). On the other hand, from Eq. (|2.4|), we have 



/^y = (Pi + Pa) = iP2 + Ph) > max{(ai + a^) , (02 + «&) }• (2.8) 

By definition, 7nx2 is the minimum of max{(ai + Oa)^, (0:2 + ttb)^} over all partitions 
of the missing transverse momentum. Therefore, we conclude that /iy > 7nx2{^N)- 
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For the reverse direction we need to prove that for a given /xat, the point {fiN, nT'T2{f^N)) 
is indeed in the consistent mass region. By the definition of mT2, there exist physical 
2+1 dimensional momenta satisfying 

2 2 2 

"l = "2 = f^N^ 

171^2 = ("1 + aaf > ("2 + 

Pl+P^2=f, Pl+Pl=f. (2.9) 

Note that if (ai + cta)^ < (0:2 + Ob)^, we can simply exchange the labels. Given ai 
and a2, we can arbitrarily choose pi^ and p2z (or equivalently, the rapidities rji, rj2) of 
particles 1 and 2, and Eqs. ( p. 91) are still satisfied. In particular, we can choose a piz 
such that r/i = rja- In this case we have {pi + PaY = (o^i + '^a)^ = my2- As for the 
other decay chain, if (0^2 + 0^,)^ = choose rj2 = rj^, if (02 + a;,)^ < have 

(P2 +Pfe)^ < when 772 = Vb, and (p2 +P6)^ — ^ 00 when 772 — ±C)0, as a result, there 
exists an ri2 such that (p2 + PbY = Tn'^2- this way we obtain physical momenta pi 
and P2 which satisfy Eqs. ( |2.7| ) with /ly = mT2{fJ'N)- This concludes our proof. 

2.3 Calculating m2-2 

In the previous subsection, we have shown that mT2 is the boundary of the mass region 
consistent with the minimal kinematic constraints. This provides us not only a way to 
understand mT2, but also an effective method of calculating it. 

We start by discussing how to determine if a mass pair (/itv, /xy) is consistent 
with the constraints in Eqs. ( |2.7|) . Note that mT2 is invariant under any independent 
longitudinal boosts for the particles a and b. This allows us to set p^ and pi to zero for 
convenience. We also assume rria > and rrib > for the moment. 

We first consider the decay chain involving particles 1 and a. From the mass shell 
constraints pf = /^^ and (pi + PaY = /^y, we can express Ei in terms of and p\: 

E. = f + ^Pl + ^y-^l'^ . (2.10) 
In order to have pi physical, we must have 

~pf = -{El - pf - pf - 1,1) < 0. (2.11) 

Eq. ( |2.11| ) imposes a constraint on possible p^ and p\. It is straightforward to show 
that the allowed {p^, p\) is the region enclosed by an ellipse. We will distinguish an 
ellipse and the region that it encloses by calling the latter an "elliptical region." The 
size of the ellipse depends on /ly monotonically. In particular, it shrinks to zero when 
f^Y = fJ'N + f^a, in which case all three particles have the same velocity. 
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Figure 2: The unbalanced solution (left) and the balanced solution (right). The red 
(dashed) ellipse and the red point are the two ellipses when /xy = fijy+nia. For the unbalanced 
solution, the point is inside the red ellipse, and mT2 = t^N + fna- For the balanced solution, 
the point is outside the red ellipse, and mT2 is given when the two ellipses (solid blue) are 
tangent to each other. 



The other decay chain is completely analogous, and we obtain another elliptical 
region for (pf, p^). However, the two decay chains are related by the measured '^t- 
Therefore, we can eliminate P^t the second elliptical region also on the 

{pf, p\) plane. In order to satisfy all the constraints, the two elliptical regions must 
overlap. Since that the two ellipses both expand as we increase /iy, mT2 will be given 
by the minimal /iy when the two elliptical regions start to overlap. To proceed we need 
to distinguish two cases, which are illustrated in Fig. ^ and discussed below. 

We assume > mb for the invariant masses of a and h without loss of generality. 
We see that we must have fJ,Y > fJ^N + f^a, otherwise the first ellipse vanishes. When 
f^Y = fJ'N + f^a, the first ellipse becomes a point, while the second ellipse has a finite 
size (or is also a point if = m^). If the point (first ellipse of zero size) is within the 
second elliptical region, then 771^2 is simply given by mT2 = f^a + f^N- This is called 
the "unbalanced configuration" in Ref. p!3| . 

The other possibility is that the point representing the zero-sized first ellipse when 
fJ'Y = fJ'N + iT^a is outside the second ellipse. In this case, we have to increase /iy until 
the two elliptical regions overlap to obtain solutions. mT2 is then given by the value 
of /iy when the two ellipses are tangent to each other. This is dubbed the "balanced 
configuration" 

Now it is clear how to calculate mT2- For a given /i^r, we first check if the two 
ellipses give us an unbalanced configuration when /ty = hn + f^a- If so, we have 
fnT2 = jJ'N + "^a- Otherwise, we need to look for the /ty when they are tangent. 
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The two ellipses are described by two quadratic equations, which can be reduced to a 
univariate quartic equation. When the ellipses are tangent, the quartic equation has 
degenerate roots and therefore its discriminant vanishes. The discriminant is in general 
a 12th order polynomial function of fi% and fiy- This provides an analytical relation 
between fi]y and 7/17-2 (/^ at) which has not be obtained in the literature before for the 
case of non-vanishing UTM. Although it would not be the most efficient way to do 
the calculation, in principle one can numerically solve the polynomial equation and 
obtain mT2- Of course, there can be more than one real solutions for the equation. 
One should keep the smallest positive /iy as 171x2 since this is the first time the two 
ellipses are tangent. 

It is convenient to use the discriminant when the UTM vanishes. In this case, the 
equations are simplified so that the 12th order equation is reduced to a 4th order one, 
for which analytical solutions are available. This confirms the existence of analytical 
solutions in the zero UTM case discussed in Ref. |T^. When UTM is nonzero, solving a 
12th order equation is numerically slow and unstable. We have developed a faster and 
more robust algorithm for calculating mT2, which is described in detail in Appendix 
0. The basic idea is that: we know that the two ellipses do not intersect when /Xy™ = 
/iAT + rria and we can also find a fiy^^ when they do intersect by an educated guess. 
Then 771^2 niust be within the interval (/iy™, /iy^^). Whether the two ellipses intersect 
can be tested easily by the Sturm sequence [1^. We repeatedly bisect the interval while 
keeping the 771^2 within it, until we reach the desired precision. 

In the above discussion, we have assumed that rria^b > 0. When either rria or rrih 
(or both) vanishes, the corresponding ellipse becomes a parabola, but the treatment 
remains the same. 



3. Mass determination using mT2 

The simplest application of mT2 is to determine my from the mT2 endpoint when rriN 
is known. However, it is often the case that m^r is also unknown, and we want to 
determine both masses simultaneously. The merit of mT2 is that it corresponds to 
the minimal constraints. Therefore it is always well-defined and calculable, which may 
prove useful at the early stage of the LHC |]T5|. If more information is available, we 
can develop more complicated methods based on mT2- 

If each decay chain in the events involves only a single two-body decay, i.e., Y 
decays to N plus a single visible particle with fixed mass, it is impossible to determine 
both masses from pure kinematics. If we consider two identical decay chains, the next 
simplest case is then that each decay chain contains two visible particles, which is 
illustrated in Fig. |^. 
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Figure 3: The event topology with 2 visible particles per decay chain. The particle X can 
be either on shell or off shell. 



When the particles in the two decay chains are identical, we can divide the situa- 
tions into four cases, depending on whether there is significant UTM, and whether the 
intermediate particle X is on shell or off shell. The latter determines whether the par- 
ticle Y decays to through a three-body decay or two consecutive two-body decays. 
For the two two-body decay cases, we would like to determine mx as well. The cases 



with zero or negligible UTM are discussed in Ref. |]10|- The authors pointed out that 
171^2^, the mT2 endpoint value as a function of fiN, has a kink at /xat = rriAr. In practice, 
it is difficult to identify the kink due to experimental smearing, but the formula for the 
m™2^ curve is known, which make it possible to fit the position of the kink. 

When the events have significant UTM, the situation is different. UTM can come 
from initial state radiation or heavier particle decays. For the latter, an example in 
MSSM is the decay chain 

q-^qX2^qfi^q^iXi, (3-1) 

where particles Xg, i and Xi are identified with Y, X and N respectively. The quark 
from squark decay can be very energetic, providing large UTM to the system. In this 
case, m™|^ curve is different from the vanishing UTM case and an analytical formula 
is in general unavailable. We focus on this case in the following, taking the process in 
( |3.1| ) as an example. 

We consider two mSUGRA points, one with > m^o and the other one with 
mi < m^o, corresponding to the three-body decay case and the two-body decay case 
respectively. The former is chosen to be the same as the model PI in Ref. |12| for 
comparison. The latter is the Snowmass SUSY point SPSla [p!6[| . 
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Figure 4: Number of consistent events as a function of (/xtvj /^y) for 1000 events, in the 
three-body decay case. The soHd Hues are the contours in 100 event intervals, beginning from 
1000 for the top contour. The dashed hne is the function /iy — iin = M_. 



1. mo = 350 GeV, mi/2 = 180 GeV, tan/3 = 20, sign(/i) = +, Aq = 0; 
m^o = 123.7 GeV, m^^ = 358.6 GeV, m^o = 70.4 GeV. 

2. mo = 100 GeV, mi/2 = 250 GeV, tan (3 = 10, sign(/i) = +, Ao = -100; 
m^fo = 181.0 GeV, m^-^ = 143.7 GeV, m^o = 100.4 GeV. 



The spectra are calculated with SPheno IjT^- We have generated 1000 events for each 
case with MadGraph/MadEvent at the parton level. In this paper, we will only consider 
the ideal case, i.e., no background, the particles are exactly on-shell, and there is no 
experimental smearing or wrong combinatorics for the visible particles. Nevertheless, 
we have avoided using features that is easily lost after the above effects are included, 
such as a kink structure, and expect our methods to be valid for realistic cases. A 
realistic study is left for a future publication. 

We first discuss the three-body decay case. Again, it is illuminating to think about 
mT2 as kinematic constraints. We have seen that for one event, the consistent mass 
region on the {^n, fiy) plane is above the curve. Using this fact, for multiple 
events we can easily count the number of events consistent with a given mass point. 
Fig. ^ is the contour plot for the number of consistent events. In particular, above 
the uppermost contour, which we identify as the m™|^ curve, the masses are consistent 
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Figure 5: Number of consistent events along the line /iy — = M_, in the three-body 
decay case. 



with all 1000 events. A kink structure in the m™2^ curve is visible. However it will be 
smeared by experimental resolutions and we avoid using it. 

Similar to Ref. |1^, we can assume that the mass difference M_ = my — rriN is 



already measured with good precision from the endpoint of the dilepton invariant mass 
distribution. Indeed, due to branching ratios, it is often the case that there are many 
more dilepton events than four-lepton events and therefore the mass difference can be 
measured much better. Drawing a line corresponding to /xy — /ijv = M_ on the (/iat, /zy) 
plane, we see that it intersects some of the contours and touches the m^f^ curve only 
at /xtv = f^N- We can draw the number of consistent events as a function of /iTv, along 
the line /xy — = M^, as shown in Fig. As expected, the number is maximized at 



fiN = mN- This is in some way equivalent to the approach in Ref. |T^, where event- 
by-event lower and upper bounds for my are obtained by intersect the mT2 curve with 
the line hy ^ fJ'N = M_. The distribution in Fig. ^ is an integral of the upper and lower 



bound distributions of Ref. |T2[ . Presenting in this way allows easy generalizations to 
other cases. To minimize the statistical error, instead of simply reading the maximum 
in Fig. ^, we can also fit the distribution to template distributions around the true uin, 



analogous to Ref. |]12 



We now turn to the two-body decay case. In this case, the particle X is on-shell, 
which gives us constraints in addition to Eqs. ( |2.7| ) (where Pa = P3+P5 and pb = p^+Pe)'- 

(Pl+P3)' = (P2+P4)' = /X^, (3.2) 

where fix is a trial mass for X. Unlike the minimal constraints, the system given by 
combining Eqs. (p^ ) and ( p.2|) can be solved event by event to yield discrete solutions 
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Figure 6: Number of consistent events distribution for = n^■T2^i^^N), in the two-body 
decay case. Left: mx is determined from the edge of the dilepton invariant mass distribution; 
right: mx is the value that maximize the number of consistent events. 



for pi and p2 (up to a four- fold ambiguity)^. Depending on whether the solutions 
are physical, we can determine if an event is consistent with a given set of masses 
(/iAT, /ix, /iy)- This fact is used in Ref. to determine all three masses by examining 
the distribution of the number of consistent events. There, the masses are obtained 
through a series of one-dimensional recursive fits. In the following, we present a sim- 
plified method utilizing the m™|^ curve. The idea is that we can follow the m™|^ 
curve, which gives us a relation between my and rriN, and count the number of con- 
sistent events. Depending on whether we want to use measurement from the dilepton 
invariant mass distribution like in the off-shell case, the method is slightly different. 

It is well known that for the decay chain in with £ on shell, there is a sharp 
edge in the dilepton invariant mass distribution at 



2, {ml - ml){ml - m%) 



'"'Uledge — 2 " V'-'-'-'/ 



^This fact was first used to study tt events in the dilepton channel, see, for example, Ref. pi 
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Assuming the edge position is measured with good precision, we obtain a relation 
among the three trial masses. Together with the relation from m™2^ curve, /ix and /ly 
are fixed for a given yU^r, up to a two-fold ambiguity from inversion of Eq. ( p.3| ). Then 
the number of events consistent with Eqs. ( p.7| ) and (|3^), as a function of //at is given 
in Fig. ^ (a). There is an evident peak at /ijy = mjy. 

The masses can be determined without using the measurement of the mu edge. 
For each fi^, we first fix /iy by the m™|^ function. Then we vary fix to maximize 
the number of consistent events. This maximum number is shown in Fig. ^ (b), as a 
function of /iiv^. Unlike the previous case, there is not a peak structure, but the number 
of events drops sharply when /iat > rriN, which can be used to estimate the masses. 

4. Discussion and conclusions 

We have demonstrated the relation between the mT2 variable and the kinematic con- 
straints for events with two identical decay chains, each of which ends up with one 
missing particle. The m™2^ curve is equivalent to the boundary of the consistent region 
in the mass space from the minimal kinematic constraints, where only the mass shell 
conditions of the decaying mother particles and the final missing particles, and the 
measured missing transverse momentum constraint are imposed. In fact, it should not 
be surprising that many different mass determination methods are closely related since 
they are based on the same kinematics. Understanding their relations may allow us 
to develop more effective and powerful ways for mass determination either by finding 
new strategies or by combining various approaches. Here we will try to give a general 
discussion of the mass determination program based solely on kinematics. 

For a given topology of new physics events, we can think of it as a map between the 
"mass" space and the "observable" space O. The mass space is the space of the mass 
parameters of the new particles which appear on shell in those events. The dimension is 
equal to the number of the unknown masses that are to be determined. The observable 
space is the multi-dimensional space of all independent kinematic observables which 
are relevant for the mass determination in those events. Basically, they are made of 
the momenta of the visible particles (jets and leptons) from the decay chains, and the 
missing transverse momentum. Therefore, each experimental event can be represented 
by a point in the observable space. In principle we can choose any basis for the observ- 
able space. It is convenient, however, to choose combinations that are invariant under 

^When m™|'^(TOjv) < my, which is always the case for finite number of events, the number of 
solvable events at /iAf = uin drops much below the maximum number. Therefore, in Fig. ^ (b) we 
have added to /xy a small constant, /iy = m™|'' + 2 GeV. 
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Figure 7: The map between a point in the mass space and the corresponding consistent 
region in the observable space. 



certain transformations which do not alter the the connection between the mass pa- 
rameters and the kinematic observables. For example, if an event of certain observable 
momenta can be produced by some mass parameters, a boost along the beam axis can 
also be produced by the same mass parameters because the momenta along the beam 
axis and the energies of the initial partons in collision are unknown. They do not have 
to be fully Lorentz-invariant as the measured missing transverse momentum breaks the 
symmetry. For two decay chains in an event which is the focus of most discussion, the 
allowed transformations are independent boosts of the two decay chains along the beam 
axis and the rotation around the beam axis, so it would be advantageous to choose the 
observable combinations that are invariant under these transformations. (If there is no 
UTM, one can also perform back-to-back equal transverse boosts on the two chains as 
the two mother particles have equal and opposite transverse momenta in this case.) 

For each point m in the mass space, there is a corresponding region /(m) C (9 in 
the observable space which is consistent with this mass point, i.e., f{m) is made of all 
possible points in the observable space that can be produced kinematically by the given 
mass parameters (see Fig. Assuming that there is a large enough set of experimental 
events from this topology with the given mass parameters at m and ignoring the issues 
such as experimental smearing and backgrounds for the moment, the region /(m) will 
be populated by these experimental events. The relative weights and densities of the 
experimental events within the region depend on other details of the underlying theory 
such as the matrix elements. On the other hand, the allowed region /(m) solely depends 
on the masses of the new particles. If /(m) is unique for each point m, then in principle 
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all the masses can be uniquely determined given enough of experimental events. It is 
possible that there are degeneracies such that different mass points map into the same 
observable region, /(m) = /(m') for m ^ m', e.g., the case of one step two-body decay 
on each chain. In that case the masses cannot be uniquely determined from kinematics 
alone and additional (model-dependent) information is required. In general we expect 
/(m) to be unique if the dimension of the observable space is large enough. From the 
above discussion, we see that the most important events for mass determination are 
those which lie near the boundary of /(m) as they determine the shape and the size of 
/(m). The edge/endpoint method can be viewed as a simple application of this idea by 
projecting /(m) down to a few one- dimensional subspaces and extract the endpoints 
of /(m) in these one- dimensional subspaces. It is also evident that it does not fully 
utilize all the relevant information contained in the experimental events as it only uses 
a few points on the boundary. In particular, in the case of two visible particles in each 
decay chain it does not give enough information to determine all masses, yet we know 
that the masses can be determined by other methods. A generalization to look at the 
boundary of the two-dimensional subspaces of /(m) is currently being studied ||l9 



It can potentially give a more powerful method than the one-dimensional endpoint 
method. Ideally, one would like to map out the whole boundary of /(m) in the high- 
dimensional observable space to get all the information contained in the experimental 
events. However, dealing with the high- dimensional space could be technically quite 
difficult. 

The method of kinematic constraints can be considered as the inverse map of the 
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mass space and the observable space discussed above. Each experimental signal event 
which is represented by a point p in the observable space O can define a region g{p) 
in the mass space which is consistent with that event (see Fig. P). Obviously, 
m G g{p) if and only if p G /(m). The correct mass point must lie within the region 
g{p) assuming that the event is a valid signal event and no experimental smearing. By 
combining many experimental events we can find the intersection of all g{pi) for all 
experimental points pi, 

G{E) = fl ^(p,), (4.1) 

Pi€E 

where E is the region in the observable space occupied by the experimental events. 
If there are sufficiently many events, we expect E = /(rritrue) for the correct mass 
point mtrue- For certain event topologies, e.g., three on-shell two-body decays for 
each decay chain |^, the intersection region shrinks to a point after combining a finite 
number of events. Obviously it would correspond to the correct mass point, G{E) = 
{rritrue}, and the mass determination is conceptually straightforward in this case. Of 
course in practice, such determination will be complicated by experimental smearing, 
backgrounds, and combinatorics. Non-trivial techniques still need to be developed 
to resolve these issues in order to demonstrate the viabihty and accuracies of mass 
determination in this way. On the other hand, for many other topologies the region in 
the mass space consistent with all experimental events G{E) remains finite. Naively one 
might think that this is a degeneracy and the masses cannot be uniquely determined. 
However, a point in the mass space m that is consistent with all events only implies 
that /(m) D E. For a generic but not the true mass point in G{E), /(m) will be 
larger than E. Assuming that the map between the two spaces is continuous, we 
then expect that its immediate neighborhood points m' in many directions will still 
be consistent with all events, /(m') D E. From this argument we see that any point 
lying in the middle of G{E) will not be the true mass point. For the true mass point, 
we should have /(rritrue) coincide exactly with E. Such a point would have the least 
neighborhood points which are still consistent with all events, or said in another way, 
it has the least degrees of freedom to move while staying within G{E). This tells us 
that the true mass point should lie on the boundary of the consistent mass region 
G{E). In particular, if there is a sharp edge or a "kink" on the boundary of G{E), 
it would be a good candidate for the true mass point. This fact has been used to 
develop a new method for mass determination in Ref. H. Now we see that the method 



of the 171x2 kink ||T0| , pJj] is another example, as the m™2^ curve is just the boundary 



of the consistent region in the two-dimensional mass (sub) space based on the minimal 
kinematic constraints. These methods effectively attempt to match the whole region of 
/(m) with the region of experimental events E as we hope to achieve in our discussion 
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in the previous paragraph. 

In the reahstic situation, such a sharp edge or "kink" of the consistent mass region 
can easily be washed out after the experimental smearing and backgrounds are taken 
into account. Therefore it may not be practical to directly search for the kink position. 
However, understanding the structure of the consistent mass region from the kinematic 
constraints allows us to develop strategies to recover the kink location by combining 
various techniques that people have developed. For example, it is well known that for 
collider signal events with missing energies, the difficulty is to determine the overall 
mass scale. The relative masses or mass differences usually can be well constrained from 
the kinematic variables such as the endpoints of invariant masses of visible particles. 
We can use those kinematic variables to reduce the mass space down to a one (or low) 
dimensional space which contains the true mass point. Then if we count the number of 
consistent events as a function of the points along this one dimensional space, it would, 
in the idealized case, exhibit a sharp turning in number of consistent events at the true 
mass point due to the "kink" nature of the consistent mass region near that point. 
Even though the sharpness of the turning point will be reduced by the experimental 
smearing and the presence of backgrounds, this "turning" feature is expected to survive 
as long as we have a reasonable data set of the signal events, and we can fit for the 
turning point to determine the overall mass scale. This was illustrated in Section ^. 

In conclusion, we have clarified the relation between the mT2 variable and the 
kinematic constraints for events with two decay chains ending with invisible particles. 
mT2 is a clever variable which simply corresponds to the boundary of the allowed mass 
region from the minimal kinematic constraints where only the constraints of mass shell 
conditions of the mother particles and the missing particles of the two decay chains, and 
the measured missing transverse momentum are used. As a by-product, we also found 
a faster algorithm to calculate mT2 from the point of view of kinematic constraints. 
These connections can also tell us how to develop new ways by combining different 
existing methods to achieve the more powerful and accurate mass determination for 
events with missing energies. It will be extremely important for reconstructing the 
underlying theory and verifying whether we have discovered the dark matter particle 
if such new physics events with missing energies are indeed found at the LHC. 
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A. The bisection method for calculating mT2 



We describe in this appendix the bisection algorithm for calculating mT2 for the bal- 
anced configuration. 

First, we need a method to quickly determine if two ellipses intersect, without 



solving the quartic equation described in Section This is done by calculating the 



Sturm sequence for the quartic polynomial, which gives us the number of real solutions 



for the quartic equation |14|. When the number of real solutions are zero, either the 



two ellipses are outside each other, or one completely contains the other one. 

For the balanced configuration, the two ellipses are outside each other for /iy™ = 
fiN + maxima, nib}. When we increase yUy, both ellipses expand. It is easy to see that 
they always intersect for /xy in some range. Thus, we need to guess a point when they 
intersect. We do this by first finding a fiy such that the two ellipses enclose a same 
point, for example, the origin. In this case, either they intersect or one contains the 
other one. If it is the former, we have found an intersecting point which is taken as 
^max_ j£ ^YiQ latter, which rarely happens, we need to do a scan from /Xy™ to find 
the intersecting point. 

After obtaining /iy™ and /iy'^^, we bisect the interval (/iy™, /iy*^^) and check if the 
two ellipses intersect at the middle point of the interval /Xy"^ = (/iy™ + Hy^^)/2. If yes, 
we set the new iiy^^ = /^y"^; otherwise, /iy™ = /iy'^. We repeat this procedure until 
the size of the interval is smaller than the precision we want. 

The algorithm has been implemented in c-|--|- and available at Ref. or from the 
authors. The code has been tested for 250k events. These include 5 datasets with 50k 
events each corresponding to the ti production in the dilepton channel, and the two 
SUSY mass points discussed in Section |^. For the SUSY points, events with UTM (from 
squark pair production and decay) and without UTM (from direct X2 P^i^ production) 
are tested separately. The tests are performed for a variety of trial masses fi^. The 
results have been compared with Ref. showing good agreement in the numerical 



values of mT2- the possibility is O(10~^ ~ 10""^) for the two programs to yield values 
that differ by 1 GeV or more, and O(10~^ ~ 10"'^) for 0.1 GeV or more. For the events 
that give small differences, our code is showing more accurate results, which can be 
verified in Mathematica by examining when the two ellipses are tangent to each other. 
Our code is also much faster (5-9 times as fast as Ref. I^H), making it advantageous 



when mT2 needs to be repeatedly calculated, for example, in evaluation of the rriTGe 



variable 13 . 
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